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Abstract. The growth of supermassive black holes (SMBH) through accretion is accom- 
panied by the release of enormous amounts of energy which can either be radiated away, 
as happens in quasars, advected into the black hole, or disposed of in kinetic form through 
| powerful jets, as is observed, for example, in radio galaxies. Here, I will present new con- 

Q j. straints on the evolution of the SMBH mass function and Eddington ratio distribution, ob- 

I i tained from a study of AGN luminosity functions aimed at accounting for both radiative 

and kinetic energy output of AGN in a systematic way. First, I discuss how a refined Soltan 
argument leads to joint constraints on the mass-weighted average spin of SMBH and of 
the total mass density of high redshift (z ~ 5) and "wandering" black holes. Then, I will 
show how to describe the "downsizing" trend observed in the AGN population in terms of 
cosmological evolution of physical quantities (black hole mass, accretion rate, radiative and 
kinetic energy output). Finally, the redshift evolution of the AGN kinetic feedback will be 
briefly discussed and compared with the radiative output of the evolving SMBH population, 
thus providing a robust physical framework for phenomenological models of AGN feedback 
£ — ' within structure formation. 
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1. Introduction: a Synthetic picture Of tion, and of the physical processes that make 

AGN evolution supernovae and gamma-ray bursts explode, the 

supermassive black holes one is the outcome 

k> I Black holes in the local universe come into two of the cosmo logical growth of structures and 

main families according to their size, as recog- of the evolution of accretion in the nuclei of 

nized by the strongly bi-modal distribution of galaxies, likely modulated by the mergers these 

the local black hole mass function (see Fig. Q). nuc i ea r black holes will experience as a result 

While the height, width and exact mass scale of the hierarchical galaxy-galaxy coalescences, 
of the stellar mass peak should be understood 

as a by-product of stellar (and binary) evolu- In the recent literature, it has become cus- 

tomary to introduce the works on cosmological 

Send offprint requests to: A. Merloni aspects of AGN astrophysics by referring to the 
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Fig. 1. The local black hole mass function, plotted as M x <f> M , in order to highlight the location and height 
of the two main peaks. The stellar mass black holes peak has been drawn assuming a log-normal distribution 
with mean mass equal to 5 solar masses, width of 0.1 dex and a normalization yielding a density of about 
1.1 x 10 7 M Q Mpc~ 3 jFukugita & Peeb les 2004) . The supermassive blac k hole peak, instead, contribute to 
an overall density of about 4.3 x 10 5 M G Mpc~ 3 jMerloni & Heinzll2008l) 



strong role they most likely play in the galaxy 
formation process throughout cosmic history. 
Indeed, a new paradigm has emerged, accord- 
ing to which the feedback energy released by 
growing supermassive black holes (i.e. AGN) 
limits the stellar mass growth of their host 
galaxies in a fundamental, generic, but yet not 
fully understood fashion. 

The strongest observational evidence for 
such a schematic picture emerged in the last 
decade. The search for the local QSO relics 
via the study of their dynamical influence on 
the surrounding stars and gas carried out since 
the laun ch of the Hubble Space Telescope 
(see e .g. iRichstone et alJU998t iFerrarese et al.1 
120081 and references therein) led ultimately 
to the discovery of tight scaling relations 
between SMBH masses and properties of the 
host galaxies' bulges (Geb hardt etalJ l2000l : 



Ferrarese & Merritt 20001: iTremaine et alj 

2002tlMarconi & Huntll2003l) . clearly pointing 



to an early co-eval stage of SMBH and galaxy 
growth. A second piece of evidence comes 
from X-ray observations of galaxy clusters, 
showing that black holes are able to deposit 
large amounts of energy into their environment 
in response to radiative losses of the cluster 
gas. From studies of the cavities, bubbles and 
weak shocks generated by the radio emitting 
jets in the intra-cluster medium (ICM) it 
appears that AGN are energetically able to 
balance radiative loss es from the ICM in the 
majority of cases (see iBirzan et al J [20081 and 
references therein). 

Nevertheless, the physics of AGN heating 
in galaxy cluster is still not well established, 
neither have the local scaling relations proved 
themselves capable to uniquely determine the 
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physical nature of the SMBH-galaxy coupling. 
As a consequence, a large number of feed- 
back models have so far been proposed which 
can reasonably well reproduce these relations. 
From the observational point of view, the cru- 
cial test for most models will be a direct com- 
parison with the high-redshift evolution of the 
scaling relations. 

There is, however, another benchmark, 
based on existing data, all models have to be 
tested upon: the evolution of the SMBH mass 
function and of the predicted energy output (ei- 
ther in radiative or kinetic form) needed to off- 
set gas cooling and star formation in galaxies. 

Here I present our recent attempt to re- 
construct the history of SMBH accretion in 
order to follow closely the evolution of the 
black hole mass function, as needed in or- 
der to test various models for SMBH cosmo- 
logical growth as well as those for the black 
hole-galaxy co-evolution. Similar to the case 
of X-ray background synthesis models, where 
accurate determinations of the XRB intensity 
and spectral shape, coupled with the resolu- 
tion of this radiation into individual sources, 
allow very sensitive tests of how the AGN lu- 
minosity and obscuration evolve with redshift, 
we have argued that accurate determinations 
of the local SMBH mass density and of the 
AGN (bolometric) luminosity functions, cou- 
pled with accretion models that specify how 
the observed AGN radiation translates into a 
black hole growth rate, allow sensitive tests of 
how the SMBH population (its mass function) 
evolves with redshift. By analogy, we have 
named this exercises AGN synthesis mod- 
elling' dMerloni & Heinzl2008l) . In performing 
it, we have taken advantage of the fact that the 
cosmological evolution of SMBH is markedly 
simpler than that of their host galaxies, as in- 
dividual black hole masses can only grow with 
time, and SMBH do not transform into some- 
thing else as they grow. Moreover, by identi- 
fying active AGN phases with phases of black 
holes growth, we can follow the evolution of 
the population by solving a simple continuity 
equation, where the mass function of SMBH at 
any given time can be used to predict that at 
any other time, provided the distribution of ac- 
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cretion rates as a function of black hole mass is 
known (see § 



In order to carry out our calculation, we 
assumed that black holes accrete in just three 
distinct physical states, or "modes": a radia- 
tively inefficient, kinetically dominated mode 
at low Eddington ratios (LK, the so-called "ra- 
dio mode" of the recent literature), and two 
modes at high Eddington ratios: a purely ra- 
diative one (radio quiet, HR), and a kinetic (ra- 
dio loud, HK) mode, with the former outnum- 
bering the latter by about a factor of 10. Such 
a classification is based on our current knowl- 
edge of state transitions in stellar mass black 
hole X-ray binaries as well as on a substantial 
body of works on scaling relations in nearby 
SMBH. It allows a relatively simple mapping 
of the observed luminosities (radio cores, X- 
ray and/or bolometric) into the physical quan- 
tities related to any growing black hole: its ac- 
cretion rate and the released kinetic power. 



In this work I will focus on just a few spe- 
cific aspects of the derived SMBH evolution, in 
particular on the redshift evolution of the mass 
function and Eddington ratio distribution, on 
the constrains we put on the mass-weighted av- 
erage spin of the SMBH population, and on the 
kinetic energy output of growing black holes. 
A more detailed discussion of the methodol- 
ogy, as well as a wi der exploration of our re - 
sults can be found in iMerloni & Heiiizl (|2008|) . 



All results will be shown accounting for the 
intrinsic uncertainties of the adopted luminos- 
ity functions. We estimated that these uncer- 
tainties can be evaluated by comparing differ- 
ent analytic parametrization of the same data 
sets; specifically, we adopted the LDDE and 
MPLE parametrizatio n for the hard X - ray lu - 
minosity function of ISilverman et all (2008), 
and two alternative parametrizations for the 
flat-spectrum radio lumi nosit y function of 
lDunlop & Peacock! dl990h and iDe Zotti et ail 
(120051) . 
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Fig. 2. Redshift evolution of the SMBH mass density (left) and average Eddington ratio (right), calculated 
for different BH mass bins (in solar mass units). In the right hand plot, black lines and grey shaded area 
represent the overall (mass-wighted) average Eddington ratio. 

knowledge of both mass, i/^M), and lumi- 
nosity, <f>({), functions, evolving the SMBH 
mass function backwards in time, up to where 
reliable estimates of the (hard X-ray se- 
lected) AGN luminosity functions are avail- 
able. The adopted hard X-ra y luminosity func- 
tion dSilverman et al.l 120081) is supplemented 
with l uminosity-depe ndent bolometric correc- 
tions ( Marconi et alj|2004l) and absorbing col- 
umn density distributions consistent with the 
X-ray background constraints (as well as the 
sources number counts, and many others), fol- 
lowin g the most rece nt XRB synthesis model 
(see lGilli etal]|2007l for details). 

In Figure [2] I show the evolution of the 
black hole mass density (left hand side) and 
of the mass-weighted average Eddington ratio, 
A = Lboi/^Edd (right hand side), both computed 
for four different black hole mass bins (in solar 
mass units). 

Between redshift zero and one, it is ev- 
ident how small mass object have a higher 
Eddington ratio, and increase their total den- 
sity much more rapidly than their high- 
mass counterpart, an effect of the well 
known phenomenon calle d "AGN downsiz- 
ing" dHeckman et al.ll2~004l) . This trend seems 
however inverted at higher redshift, when the 



2. The evolution of SMBH Eddington 
ratio 

We studied the evolution of SMBH mass func- 
tion through a continuity equation that can be 
written as: 

d J^ + ±( H M,t)(M(M,t») = (1) 

where M is the black hole mass, /u = Log M, 
jj. = Log M, i//(M, t) is the SMBH mass func- 
tion at time f, and (M(M, t)) is the average ac- 
cretion rate of SMBH of mass M at time f, 
and can be defined through a "fueling" func- 
tion, F(fi, fi, t), describing the distribution of 
accretion rates for objects of mass M at time 
t: (M(M,z)) = J MF(Ji,fi,z)dfi. 

Such a fueling function is not a priori 
known, but can be derived inverting the inte- 
gral equation that relates the luminosity func- 
tion of the population in question with its mass 
function: 



cf>(£,t) = J F(U,n,t)ittn,i)&n (2) 

where I have called £ = LogL\, Q \. 

Thus, we have integrated eq ([TJ starting 
from z = 0, where we have simultaneous 
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largest black holes are assembled; a more pre- 
cise assessment of this phenomenon will re- 
quire better statistics on the high redshift AGN 
luminosity functions. 

3. Constraints on the average black 
hole spin 

One of the most far-reaching conclusions of 
the analysis presented in Merloni and Heinz 
(2008) was the necessity of a broad distribu- 
tion of Eddington ratio for SMBH of vari- 
ous masses and redshift, in order to reconcile 
the observed evolution of the AGN luminosity 
functions with that of the mass function itself. 
Consequently, there are always black holes of 
any mass accreting in different "modes" at any 
time, and this bears some consequences on the 
final estimates of the average efficiencies. 

In fact, in our three-mode scheme for black 
hole accretion, the radiative efficiency, e ra d, 
is not necessarily equal to the accretion effi- 
ciency, 77. The latter represents the maximal 
amount of potential energy that can be ex- 
tracted, per unit rest mass energy, from mat- 
ter accreting onto the black hole. This quantity, 
77(a), depends on the inner boundary condition 
of the accretion flow only, and, in the classical 
no-torque inner boundary, is a function of BH 
spin, a, only, ranging from rj(a — 0) - 0.057 
for Schwarzschild (non-spinning) black holes 
to T](a = 1) st 0.42 for maximally rotating Kerr 
black holes. 

On the other hand, the radiative efficiency, 
frad = Lb \/Mc 2 depends both on the accretion 
efficiency (i.e. on the inner boundary condi- 
tion of the accretion flow) and on the nature of 
the accretion flow itself. Based on our current 
knowledge of the physical properties of low 
and high lumin osity AGN (and s t ellar m ass 
black holes), in iMerloni & Heinzl d2008l) we 
have adopted the following parametrization: 

r, ■ \ I 1. m>m a ... 

frad = TlJ(tn) = J1X<,.,, » . . . (3) 

'■' v ' 1 I (m/m cr ), m < m cr 

where m cr = A cr is the critical Eddington- 
scaled accretion rate above which the disc be- 
comes radiatively efficient, assumed, in our 
computation, to be at the universal value of 
0.03. 



In lMerloni & Heinzl d2008l) we have shown 
how using standard I Soltanl dl982h type of ar- 
guments, i.e. comparing the local mass density 
to the integrated mass growth in AGN phases, 
very tight constraints can be put on the average 
radiative efficiency of the accretion process: 



ft(i + 6ost) ~ 6(1-6 + 6»t)' 

where £o is the local mass density in units of 
4.3 xl0 5 M o Mpc~ 3 , while £ and £ os , are the 
mass density of z ~ 5 and "wandering" SMBH, 
respectively (also in units of the local mass 
density). 

In order to translate our constrains on the 
radiative efficiency onto a constrain on the av- 
erage black hole spin, however, we need to 
carefully consider the distribution of angular 
momentum of the matter accreting onto th e 
black holes. As discussed in Kin g et al.1 (2008), 
the most important factor is the amount of mass 
that can be accreted "coherently", i.e. keep- 
ing the same large scale angular momentum. 
If this is larger than a few per cent of the 
black hole's mass, the black hole will be al- 
ways spun up, and most of the SMBH mass 
is accreted from a prograde disk. If, on the 
other hand, accretion proceeds via small, in- 
dependent (randomly oriented) sub-units, pro- 
grade and retrograde a ccretion events a re al- 
most equally probable dKing et al.ll2008l) . and 
we should use a "symmetrized" relation be- 
tween black hole spin and accretion efficiency: 
T] ± (a) = (l/2)[T)(a) + T)(-a)]. Most likely, the 
true average relation between spin and effi- 
ciency will lie somewhere in between these 
two extreme cases, depending on the detail of 
the SMBH fueling mechanism at various red- 
shifts. 

Figure [3] shows the constraints derived in 
the ^i ost , (a*) plane, where (a*) is the mass 
weighted average spin parameter of the SMBH 
calculated in verting e i ther t he classical GR 
77(a) relation d Shapiro! (120051) : left panel, pro- 
grade accretion only) and its symmetrized ver- 
sion (right panel, random accretion), and tak- 
ing into account the assumed relationship of 
eq. ([3j between 77 and 6^. 

It is interesting to note that both the amount 
(numbers and masses) of black holes effec- 
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Fig. 3. Constraints on the mass-weighted average spin of SMBH and on the fraction of "wandering" black 
hole mass density £ tost . The dark blue contours between solid lines correspond to the case of negligible black 
hole mass density at z = 5 = 0), progressively lighter blue areas represent the case of larger and larger 
values of The left hand side shows the calculation assuming accretion onto SMBH proceed always in 
prograde fashion; the right hand side shows the opposite extreme of purely random accretion. These two 
extreme cases encompass all possible cosmological solutions. 



tively ejected from galactic nuclei due to grav- 
itational wave recoil after a merger and the 
average radiative efficiency of accreting black 
holes depend on the spin distribution of evolv - 
ing SMBH (see e.g. iBerti & Volonterill2008l) . 
Thus, the plots of Fig. [3] couple implicitly 
the spin distribution of accreting black holes 
(through fi ost and e la d) and the properties of 
the seed black hole population, whose density 
must be reflected in the z ~ 5 mass density 

4. The kinetic energy density output 
of SMBH 

Finally, I discuss here briefly the consequences 
of our synthetic picture of AGN evolution for 
the issue of AGN feedback in the form of ki- 
netic energy associated to radio jets. 

In order to do that, we start from the 
recently found correlation between kinetic 
power and radio core lumin osity of AGN 
jets dMerloni & Heinzl 120071) . and use the 
flat sp ectrum radio luminosity f unction evolu- 
tion dDunlop & Peacockl Il990t |Pe Zotti et al.l 
120051) to derive the evolution of the kinetic 



luminosity function of AGN. The effects of 
beaming have been taken into account statisti- 
cally, bofh_m_fh£_ccre-kinetic luminosity rela- 
tion dMerloni & Heinzll2007l) and in the evolv- 
ing flat spectrum radio luminosity functions 
dMerloni & Heinzll2008l) . 

By integrating over the kinetic luminosity 
functions, we get an estimate of the AGN ki- 
netic energy density output as a function of 
redshift, as shown in the left hand panel of 
Fig. 21 where a comparison is made with the 
(bolometric) radiation energy density output 
and the radiation energy density in the radio 
band (5 GHz). The most notable feature of 
this plot is the markedly different redshift evo- 
lution of the radiative and kinetic power out- 
put, with the latter showing a much smaller 
amount of evolution between z — and z — 1 . 
This is due essentially to the increasing num- 
ber of SMBH accreting at low Eddington ra- 
tios in the so-called "low-kinetic" (LK) mode, 
where a substantial fraction of the gravitational 
energy of the accreted mass is converted into 
jet kinetic power. Such a wea k redshif t depen- 
dence, already suggested by iMerlonil (120041) . 
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is in fact necessary in order to reproduce 
the high-end of the galaxy mass function in 
semi-analytic models of s tructure formation 
that invoke AGN feedback dCroton et alJ200d 
iBower et al.ll2006l) . 

Let us now examine the kinetic energy 
production efficiency of growing black holes. 
We first compute the total integrated (mass 
weighted) average kinetic efficiency as^: 

fr Pkm(z)dz 

<e k i„> = 7§- (2.8 ± 0.8) x 1(T 3 (5) 

where ^P^iz) is the total black hole accretion 
rate density at redshift z. Combining the re- 
sults discussed in section [2] on the average ra- 
diative efficiency with those shown in eq. © 
we conclude that SMBH, during their growth 
from z ~ 5 till now, convert about 25 (down to 
15 if hidden cores are considered) times more 
rest mass energy into radiation than into ki- 
netic power, with the exact number depending 
on the poorly known details of the intrinsic jet 
cores luminosity function, as well as on our as- 
sumptions about the beaming corrections to be 
made. Similar, but complementary studies of 
the AGN kinetic luminosity function evolution 
based instead on steep spectrum sources ha s 
been recently c arried out lKording et al.l(l2008l) : 
IShankar et al.l (|2008|) . Reassuringly, the results 
of both these works are in reasonable agree- 
ment with those I presented here. 

We can also compute directly the kinetic 
efficiency as a function of redshift and SMBH 
mass today, which I show in the right hand 
panel of Fig. |U with the horizontal solid 
lines showing the mass weighted average from 
eq. (O. The various curves describe the main 
properties of kinetic feedback as we observe it. 
For each of the chosen mass ranges, the kinetic 
efficiency has a minimum when black holes of 
that mass experience their fastest growth: this 



1 For the sake of simplicity, I discuss here only 
the case in which we derive the intrinsic radio core 
luminosity function from the flat spectrum sources 
only; an upper limit to the contribution from hid- 
den cores in steep spectrum radio sources have been 
considered in Merloni and Heinz (2008), yielding a 
final upper limit for the kinetic efficiency about a 
factor of 2 higher 



is a different way of restating the conclusion 
that most of the growth of a SMBH happens 
during radiatively efficient phases of accretion. 
However, when the mass increases, SMBH are 
more and more likely to enter the LK mode, 
which increases their kinetic efficiency. More 
massive holes today, have entered this phase 
earlier, and by z = they have reached the 
highest kinetic efficiency of a few xlCF 2 . This 
is a natural consequence of the observed anti- 
hierarchical growth of the SMBH population, 
and of the chosen physical model for the ac- 
cretion mode of low-Eddington ratio objects. 

5. Conclusions 

I have outlined some recent results of our work 
aimed at pinning down as accurately as possi- 
ble the cosmological evolution of active galac- 
tic nuclei and of the associated growth of the 
supermassive black holes population. 

In particular, I have focused here on the 
global (integrated) constraints on the mass- 
weighted average spin of SMBH, and I have 
discussed in some details the specific ways in 
which these are tighten to the very interesting 
open issues regarding the population of high- 
redshift SMBH and that of black holes ejected 
from galactic nuclei due to gravitational wave 
recoil in merger events. 

I have also discussed a few generic prop- 
erties of the kinetic energy output of growing 
black holes, emphasizing the importance of a 
late phase of low-luminosity, jet-dominated ac- 
cretion onto the most massive objects. 

The richness of details we have been able 
to unveil demonstrates that times are ripe 
for comprehensive unified approaches to the 
multi-wavelength AGN phenomenology. At 
the same time, our results should serve as a 
stimulus for semi-analytic and numerical mod- 
elers of structure formation in the Universe to 
consider more detailed physical models for the 
evolution of the black hole population. 
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Fig. 4. Left: the AGN kinetic energy density output as a function of redshift, (light grey area between 
dashed lines) compared to the (bolometric) radiation energy density output (grey area between solid lines) 
and the radiation energy density in the radio band (5 GHz; dark grey area between dot-dashed lines). Right: 
Redshift evolution of the kinetic efficiency. SMBH of different masses at z = are here plotted separately, 
with a color coding analogous to that of Fig. [2] the horizontal black solid lines mark the mass weighted 
average values for the kinetic efficiencies, with the dashed lines representing the uncertainties from the 
particular choice of radio and X-ray luminosity functions (see Merloni and Heinz 2008 for more details). 
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